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. The radiative Higgs decays h — )• ^I'^l^ with I = e, fj, and r are analyzed in the stan- 

. dard model using mh = 125 GeV. Both tree and one-loop diagrams for the processes 

are evaluated. In addition to their decay rates and dilepton invariant mass distribu- 
tions, we focus on the forward-back asymmetries in these modes. Our calculation 
^ . shows that the forward-backward asymmetries in /i — )• je'^e~ and h — )■ ji-i~^fi~ could 

■ be up to 10~^ while in the t~^t~ final state, these asymmetries are below 1%. Thus 

the forward-backward asymmetries in /i — t- ^l~^l~ might be interesting observables 
in the future precise experiments both to test our understanding of Higgs physics in 
the standard model and to probe the novel Higgs dynamics in new physics scenarios. 
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1 Introduction 



Since the mass term in the gauge theory violates the gauge invariance, gauge bosons 
in SU{2)l X U{1)y gauge field theory should be massless. The Higgs mechanism [1] 
provides a quite simple approach to circumvent this requirement. One can begin with 
a gauge invariance theory, which has massless gauge bosons and the unstable vacuum, 
and obtain a theory with massive gauge bosons and the stable vacuum after the algebraic 
transformation on the Lagrangian due to the spontaneous symmetry breaking of the 
vacuum. The simplest form of this mechanism realized in the standard model (SM) 
requires the existence of a single neutral Higgs boson. Since all masses of massive particles 
in the SM are originated from their couplings with Higgs and the Higgs sector plays a key 
role in our understanding of the nature of the world, the searching for the Higgs boson 
was one of the most important tasks in the past few decades, and was also one of the 
main motivations for the construction of the Large Hadron Collider (LHC) at CERN. 
After more than forty years' efforts, we finally found a Higgs-like scalar particle at around 
125 GeV recently thanks to the hard work of ATLAS [2] and CMS [3] Collaborations at 
LHC. 

Theoretically, the property of the SM Higgs boson has been extensively studied. The 
experimental analysis of several different decay channels of the new particle, observed by 
ATLAS [2] and CMS [3] , shows that its properties are consistent with the elementary Higgs 
boson in the SM. However, some unexpected signals also appear. The excess events of 
the diphoton channel /i — )■ 77 have drawn many attentions, which have been investigated 
in a large number of papers using various models [H O [6] . Due to the limitation of the 
present experimental data, it is still a long way to distinguish these different scenarios. 

On the other hand, besides the /;.—)■ 77 decay, a complementary channel, the radiative 
Higgs decays h — )■ {I = e, fi and r) could also provide some useful information on 

the SM Higgs boson. Actually, these radiative decays have been studied in the literature 
[3 El El Uni [H] and the decay rates and invariant mass distributions have been evaluated. 
In the last two papers [101 E], the newly reported mass value of the Higgs boson candidate 
from pi E] has been used in their calculations. The main purpose of the present paper is to 
focus on the angular distributions, in addition to the rates and invariant mass distributions 
of the decays. It is known that, in the SM, h — >■ 7/+/" decays receive the different type 
contributions from both tree and loop diagrams [H El [10] . The angular distribution, due 
to the interference of these different contributions to the decay amplitudes, will lead to a 
quantity, called as the forward-backward (FB) asymmetry, which might be an interesting 
observable in the future precise experiments to further understand the properties of SM 
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Figure 1: The tree diagrams for h — ^yl^l 



Higgs Boson, as well as to explore new physics in the Higgs sector [121 1T3] . 

The paper is organized as follows. In section 2, the amplitudes of the process h — 
'-yl^l' from the different diagrams are evaluated, and the expressions of the differential 
decay rate and the FB asymmetries are presented. Section 3 is our numerical analysis and 
some discussions on phenomenology. Finally, we summarize our conclusions in section 4, 
and explicit expressions of some loop functions are given in the Appendix. 



2 The amplitude oi h ^ 

The tree diagrams of the processes h are not forbidden, where the photon can 

be radiated from the lepton leg after Higgs decay to /*/, as shown in Figure [U Due to 
the helicity suppression, this contribution is proportional to the mass of lepton. It will be 
shown that the contribution of the tree diagrams can be neglected safely in the process 
h — 7- 7e+e~ due to the large mass hierarchy between electron and Higgs, but for ^~ 
and '-/r^r^ in the final states, the tree diagrams should be included. This is consistent 
with the results in [71 El HQ] . 

From Figure [H the decay amplitude of /i — t- at the tree-level can be expressed 

as follows 

where Cq = — J^^II^q^ , f^^i is the mass of lepton, is the fine-structure constant, Ow is 
the electroweak mixing angle, and /ci, ^2 and p represent the momentum of /~ and 7 
in the final states, respectively. 



2 



w 



7 



7 



7 



Figure 2: The photon and Z pole one-loop diagrams for h — t- 7/+/ . 
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Figure 3: Four-point box diagrams for h 'yl'^l . The photon can also be radiated from 
another W line in the left diagram and from other charged lepton lines in both diagrams. 

Since its tree-level amplitude is proportional to the mass of lepton, the next-to-leading 
contributions could be very important for the h — '-fl~^l~ transitions, especially for the 
electron and muon modes. The typical one-loop Feynman diagrams for these processes 
have been displayed in Figures [2] and [31 respectively, which are of two basic types: (i) 
the photon and Z pole diagrams via h — 77* — )■ '~fl~^l~ and h — '-/Z* — )■ '~fl^l~ (Figure 
|2]); and (ii) four-point box diagrams involving virtual W and Z gauge bosons inside the 
loop (Figure [3]). It has been pointed out by the authors of Ref. [7] that type (i) gives the 
dominant contributions. Our calculations give the consistent conclusion with [7]. One 
can check Figures HI [5l and [6] in the next section for the numerical analysis, that type 
(ii) diagrams have only very little contributions. In the remainder of this section, for 
simplicity, we only show explicitly the results from Figure [2l to illustrate the discussions 
on the differential decay rates and the FB asymmetries. Actually, we include the four- 
point box diagrams' contributions in our numerical calculations of section 3. 
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The amplitude of h -fl^l at the one-loop level can be expressed as 
^ML = te''*[ip,qu-9f.uP-q)u{k2){Cir + C2ri'Hh) 



where 



1 . o . \ „ „ 1 „ ^ 1 



Ci — — ( ^ ~ sin 9w j Pz^s-yz — ^n^7 5 ^2 — -Pz^s^Z, 

Cs = - ( ^ - sin^ j PzUa^Z, Ci = ^Pzlla^z (3) 



with 



1 1 



(4) 



sin 6w cos 6w q"^ — fn\ + '^'^z^z 

Here we denote q as the momentum of the virtual particle such as 7* or Z* in Figure [21 
and = {ki + ^2)^ is lepton pair mass squared. Ha-yz, ^s-yz and 11^^, induced from the 
Z* and 7* pole diagrams in Figure [21 respectively, are given by 

Ha^z = ^ -an -4/2 (r/, A/), (5) 

niw sm Bw sm tlw cos 6w 

al 

Hs-yz = ^-7; — [— cot 9wAw {tw, Xw) 

mw sm tlw 



sm fe^vy COS 



(6) 



2 

n^7 = [-^w {rw, Xw) - 4:N,Q}Afi (r/, A/)] . (7) 

777,14/ sm U\Y 

Here Aw, Afi, and Af2 are functions denoting the contributions from different loops (W 
loops and fermion loops), with Tj = Aj = ^^(i = /, VT), and is the mass, is 
the color multiplicity, Qf, in unit of e, is the charge, Tf is the third component of weak 
isospin of the fermion / inside the loop. Expressions for loop functions A^'s have been 
shown in the Appendix. 

It is easy to see that the 7* and Z* pole loop diagrams of Figure [21 are direct extension 
of the processes /i — 77 [Hj and h — j- 7Z [15]. One can recover those results of [TH 
\TE\ by setting 7* and Z* on-shell in the above expressions excluding the Ha^z parts. 
The Ha-yz part amplitudes, proportional to Af2, receive contributions from the fermion 
loop diagrams in which Z gauge boson is axial-coupling to the loop fermion. This will 
generate a parity-odd h — )■ 7Z amplitude. As is well-known that the h — ?■ 7Z transition 
is dominated by W^-loop contribution, which only leads to the parity-even amplitude. 
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The small parity-odd amplitude can be safely neglected in the calculation of r(/i — )■ 7Z) 
since it does not interfere with the dominant parity-even amplitude. However, it will be 
a different story for the h — )■ '-fl~^l~ decays because the dominant parity-even h'-fZ vertex 
can also generate the parity-odd h — )■ '~fl^l~ amplitude through the axial-coupling of Z 
gauge boson and leptons in the final states. This means that this parity-odd amplitude 
cannot be simply abandoned in the present work. Such structures have also been shown 
in Ref. [10]. More interestingly, it will be shown below that, these structures {Ha^yZ part 
amplitudes) play a key role in the FB asymmetries of the h — j- decays in the SM. 

The differential decay rate of /i — )■ '~fl^l~, including both tree and loop diagrams 
contributions, can be written as 



dV _ {ml - 
dq'^dcosO 5127r'^m^ 



A[|Co|M + 2Re{CoC*)B + 2Re{CoCl)C 



h 

+ \Cs\^)D + + \C,\^)E + 2Re{C,C2 + CaC*)^] (8) 



with 



16 

ml - q^y(r- cos ey 
-{mi + q'-8mlq')Pfcos'e], (9) 



A = ^ _ o^2.7_ 02 n^2 K + + 32mf - Sm^q^ - Sm^ml 



D = ^'^'-^'^\ q^ + 4m^ + q^/3fcos'9), (12) 

E = ^^^i^g^A^(l + cos^^), (13) 
F = {ml-qyq^f3iCOs9, (14) 



where A = y 1 ~ ^ is the angle between the three momentum of Higgs boson and 
the three momentum of /~ in the dilepton rest frame, and the phase space is given by 

Amf <q^< ml, -1 < cos6' < 1. 

The contribution of the tree diagrams is given by A term in eq. (|8]) and this term will 
be divergent when q^ — > m^, which corresponds to the soft photon in the final states. 
To avoid such divergence, we cut the minimal energy of photon at 1 GeV when we carry 
out the integration over q^. B and C terms show the interference between the tree and 
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one-loop diagrams and they are both suppressed by the mass of lepton. The last three 
terms E, and F, originate from the contribution of the one-loop diagrams. 

Although C and F terms, which are proportional to cos^, do not contribute to the 
decay rate, these two terms will lead to a very interesting observable, called as the FB 
asymmetry, which is defined as 

:dcos6 — / — ^ dcos6 



A / ■z\ dq^d COS 9 J_i da d cos 9 

^FB(g ) = — — To Jf. , (15) 

dcos 9 + — ^ acos 9 



dq dcos9 J^i dq dcos9 

and the corresponding integrated FB asymmetry over is 



9max ri? 



dT /'^■"^'' f dT 

— 2 dcos9dq^ — / / — 5 dcos9dq^ 

Ja'^ Jo dq dcos9 Ja^ i_i dq dcos9 , , 

-^FB - —2 -j — -2 — ) IJ-Oj 



dcos9dq^ + / / — 5 dcos9dq 



dq dcos9 Jq2, J_i dq dcos9 



2 

min ' " ^ ^min 



2 



where {q^^i^, 5'max) denotes the range of the integration over g^. From eq. (|H]), one can thus 
find that the nonzero C3 or C4, which is proportional to Tla^yz, is required to generate the 
nonzero FB asymmetries in the SM. Since C term is proportional to mi, it is expected 
that F term would give the dominant contributions to the FB asymmetries in /i — )■ 7e^e~ 
and h — > 7/i"*"yU^ decays. 

We adopt different kinematical variables for the differential decay rate from those in 
Refs. [TOj [TT] because we are more concerned about the FB asymmetries in the present 
work. As we shall see, the cos^ dependence in the differential decay rate could provide 
some interesting information in h ^ decays. 



3 Numerical analysis 

3.1 Decay rates and invariant mass distributions of /i — > 7/+/ 

Both tree and loop diagrams can contribute to the h — > decays, and the amplitudes 

at the tree-level are proportional to the mass of the lepton in the final states. Our 
numerical calculation shows that, for m/j = 125GeV in the SM, 

T{h ^77) 
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Figure 4: The invariant mass distributions of /i — ?■ 7e"'"e~ normalized by r{h — i- 77). 
The red hne denotes the contribution of the tree diagrams, the thin sohd hne denotes the 
contribution from the 7* pole diagrams, and the dashed line the contribution from the Z* 
pole diagrams while the thick line gives the total contributions. The dotted line denotes 
the contribution from the four-point box diagrams. 



for the electron mode, and the contribution of tree diagrams is vanishingly small; for the 
muon mode, 

and its tree diagram could play a relevant role, which gives about 30% contribution to 
the rate; while in the process h — )■ '-fT~^T~ , we have 

— = 3.04, 19 

r(/i^77) ^ ^ 

where the tree-level contribution is dominant, and the loop diagrams give about 1% 
contribution. As mentioned above, in order to avoid the infrared divergence from the tree 
diagrams, we cut the minimal energy of photon at 1 GeV to illustrate our results of f|T7|) . 
(HHD and (USD. 

The dilepton invariant mass distributions oi h ^ normalized by r{h 77), 

have been given in Figures HI [5l and El respectively. Different types of contributions in- 
cluding the tree diagrams, 'y*/Z* pole diagrams, and four-point box diagrams, are plotted 
separately for comparison. It is shown that the four-point box diagrams indeed give a 
negligible contributions for all three modes, and 7*/Z* poles (corresponding to the two 
peaks) give the dominant contributions to the invariant mass distribution in h ^ '-fl^l~ 
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Figure 5: The invariant mass distribution of /i — )■ ^fi^fi~ normalized by r(/i — )■ 77). The 
red hne denotes the contribution of the tree diagrams, the thin sohd hne denotes the 
contribution from the 7* pole diagrams, and the dashed line the contribution from the Z* 
pole diagrams while the thick line gives the total contributions. The dotted line denotes 
the contribution from the four-point box diagrams. 
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Figure 6: The invariant mass distribution of /i — )■ 7r^r~ normalized by r{h — 77). The 
red line denotes the contribution of the tree diagrams, the thin solid line denotes the 
contribution from the 7* pole diagrams, and the dashed line the contribution from the Z* 
pole diagrams while the thick line gives the total contributions. The dotted line denotes 
the contribution from the four-point box diagrams. 
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Figure 7: The FB asymmetry Ap^lq"^) in /i — )■ 'ye~^e~ as a function of g^. The dotted 
hne denotes the contribution from the 'j*-Z* interference, and the dashed hue denotes the 
contribution from the Z*-Z* interference while the sohd hne gives the total contribution. 
The contribution through the interference between the tree diagram and Z* pole diagrams 
is extremely small to be shown here. 

with I = e, /i. From these three plots, one can readily understand our above results of 
the decay rates, that the tree diagram contribution is vanishingly small in the e'^e~ case, 
relevant in the fi~^fi~ case, and dominant in r+r" channel. Note that these invariant mass 
distributions in /i — > •ye'^e' and h — )■ 7/i"'"/i~ have been studied in Ref. [H] and Ref. p!0] . 
respectively. One can find that our plots (Figures 6 and 7) are basically consistent with 
those results. However, because of the different cuts and different normalization, it is not 
easy to perform a detailed comparison. 

3.2 Forward- backward asymmetries in h ^ 

As expected in section 2, the FB asymmetries ^fb(q'^) in /i — > 7e+e~ and h — '-ffi~^fi~ 
decays mainly arise from the F term in eq. ([8]), which, by recalling eq. ([3]), consists 
of the contributions from the interference between the 7* pole amplitude and the Z* 
pole amplitude (referred as 7*-Z* interference), and from the interference among the 
different amplitudes of Z* pole diagrams (referred as Z*-Z* interference); while C term 
contribution generated from the interference between tree and loop diagrams could be 
neglected or very small. This is confirmed by our numerical calculations. 

The FB asymmetries Ap^iQ"^) in /i — )■ 7e+e~ and h — )■ 7/1 """/i" have been plotted in 
Figures [7] and m respectively, which, in most kinematical region, could be up to 10~^. One 
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Figure 8: The FB asymmetry Afb(q'^) in /i — ?■ '-ffi'^fi~ as a function of g^. The dotted 
hne denotes the contribution from the 7*-Z* interference, and the dashed hne denotes 
the contribution from the Z*-Z* interference, and the red hne denotes the very small 
contribution through the interference between the tree diagram and Z* pole diagrams 
while the solid line gives the total contribution. 
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Figure 9: The FB asymmetry AfbIq'^) in /i — t- 7r+r~ as a function of g^. The dotted line 
denotes the contribution from the 7*-Z* interference, the dashed line denotes the contri- 
bution from the Z*-Z* interference, and the red line denotes the contribution through the 
interference between the tree diagram and Z* pole diagrams while the solid line gives the 
total contribution. 
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102-302 


302-502 


502-702 


702-902 


902-1102 


full phase space 


h — 7e"'"e~ 


-0.5% 


-2.3% 


-3.7% 


0.8% 


1.9% 


0.4% 


h — )■ 7/x"'"yLi~ 


-0.4% 


-1.9% 


-2.3% 


0.7% 


1.5% 


0.4% 



Table 1: The integrated FB asymmetries ^fb for some cuts on in h ^ 'ye'^e and 
h — )■ 'jfJ^'^fi^ decays. 



can see that the C term contribution is very small for the muon mode, and we even do 
not plot this part contribution in Figure [7] since it is extremely small to be shown. Figure 
Ogives the FB asymmetry in /i — t- ^t^t~ , which is about one order of magnitude smaller 
than the former two modes. It is thought that C term contribution in eq. ([8]) might be 
important for r channel, however, the tree-level amplitude is strongly dominant over the 
other amplitudes, which, from eq. (IT^ . obviously leads to the suppression of ApB{q'^) in 
h — !■ jT^r^ . 

As shown in Figures [7] and [H Ap^lq"^) in /i — > 7e"^e~ and h — ?■ '-ffi~^fi~ will be from 
negative to positive when g2 jg running in the phase space. This comes from that the 
7*-Z* interference will change sign for q"^ crossing m|, since, in this case, the real part 
of Pz in eq. (jlj) changes sign. Thus the integrated FB asymmetries ^fb defined in eq. 
(fT6|) may be not very significant if we integrate q^ in the full phase space. On the other 
hand, when we calculate these integrated asymmetries for some cuts on q"^, Afb could be 
at percent level. The corresponding numerical results have been given in Table [H and 
^FB in the full phase space of q"^ is only 0.4% for both modes. While for some kinematical 
regions, these asymmetries might be enhanced, for instance, 50 GeV < < 70 GeV, 
the magnitude of ^fb is 3.7% for h — •ye'^e' and is 2.3% for h — 7/i+yU~. In the process 
h — )■ '~fT~^T~, the integrated FB asymmetry will be also very small, .4.fb~ 0.02% in the 
full phase space of g2. Therefore it is not significant to calculate them for any cuts on q'^. 

4 Summary and Conclusions 

We have presented a comprehensive analysis of the radiative Higgs decays h '~fl^l~ 
with / = e, /i and r in the SM, by using mh=l25 GeV, the newly reported mass value of 
the Higgs boson candidate from LHC [21 13] . Both tree and one- loop Feynman diagrams 
for the processes are evaluated. It is found that loop amplitudes are dominant for both 
electron and muon modes. The tree-level contribution can be fully neglected in /?, — )■ 7e"*"e~ 
decay, and for the muon mode, it will give a sizable contribution to the decay rate while 
the tree-level transition is very important in h '~fT~^T~. Our numerical results show 
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that T{h — jl'^l^) with I = e, is about 6% of T{h — )• 77), and T{h — > 7r+r") is 
about factor three of the diphoton rate, which is almost entirely contributed by the tree 
diagrams. The dilepton invariant mass distributions, normalized by r(/i — 77), have 
been analyzed for all three modes. Different types of contributions are plotted separately 
for comparison, the four-point box diagrams are found to be negligible, and the photon 
and Z pole diagrams give the dominant contribution in /i — 7e+e~ and h — '-ffi'^fi'. Our 
results are basically consistent with those of Ref . [10] and Ref . [H] , in which the different 
cuts and normalization have been used. 

One of the main motivation of the present work is to investigate the FB asymmetries 
in h ^ Different from the two-body decay h 77, three-body radiative decay 

may generate a nontrivial angular distribution, which would provide some complementary 
information for the diphoton decay. In the SM, thanks to the parity-odd h — > 7Z* 
amplitude induced from the fermion loop, nonzero FB asymmetries can be expected, and 
we have explicitly evaluated the SM contributions to Ap^iQ"^) and ^fb in these decays. 
Interestingly, these asymmetries can be up to 10~^ in /i — > '~fe~^e~ and h — '~ffi~^fi~ decays 
while they are suppressed below 1% for the 7r^r~ final states. Consider the magnitude 
of these asymmetries, the measurement of them is indeed not an easy task; nevertheless, 
experimental studies of the FB asymmetries in /i — j- 7/"*"/" decays in the future would be 
very helpful both to increase our understanding of the properties of SM Higgs boson and 
to probe the novel couplings of Higgs in the new physics scenarios. 

In summary, our present analysis indicates that the radiative Higgs decays h — > '-yl^l~ 
are worth to be seriously considered in the future experimental investigations. These 
decays, due to their complementarity to the /i — > 77 decay, may provide us with some 
valuable information on the Higgs physics in the SM and its possible extensions. 
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Appendix: Loop functions 



Here we explicitly show the expressions for some loop functions in section 2. 
AJ^{r,\) = /i(r,A)-/2(r,A) , 
^/2(r,A) 



A 



[2g{r) - 2g{\) + f{r) - /(A)] 



with 



/i(r,A) 
/2(r,A) 



where 



Aw{t, A) 



-16(l--)/2(r,A) 



/i(r,A) 



rA 



+ ./'^'.. [/(r) - /(A)] + T^^W) - g{\)] 



2(r-A) 2(r-A)2 
rA 



(r - A)^ 



2(r - A) 
f{r) = 
9{r) = 



[fir) - /(A)] 



arcsm 



2 1 



1 2 



Jr — \ arcsin ^ 



In 



r > 1 
r<l , 

r > 1 
r < 1. 



(A.l) 
(A.2) 



(A.3) 

(A.4) 
(A.5) 

(A.6) 



(A.7) 



For ruh = 125 GeV, r < 1 will not be used if we only consider W and top quark loop. 
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